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The shear viscosity of an interacting pion gas is studied using the Kubo formalism as a microscopic 
description of thermal systems close to global equilibrium. We implement the skeleton expansion in 
order to approximate the retarded correlator of the viscous part of the energy-momentum tensor. 
After exploring this in g<f> 4 theory we show how the skeleton expansion can be consistently applied 
to pions in chiral perturbation theory. The shear viscosity r\ is determined by the spectral width, 
or equivalently, the mean free path of pions in the heat bath. We derive a new analytical result for 
the mean free path which is well-conditioned for numerical evaluation and discuss the temperature 
and pion-mass dependence of the mean free path and the shear viscosity. The ratio r\js of the 
interacting pion gas exceeds the lower bound 1/4-7T from AdS/CFT correspondence. 
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I. INTRODUCTION 

One of the lessons learned from experiments at the 
Relativistic Heavy Ion Collider (RHIC) in search for the 
quark-gluon plasma is the fact that the deconfined quark- 
gluon matter behaves as an almost-perfect fluid HHZ]. 
Nowadays we can look forward to data from the Large 
Hadron Collider (LHC) at even higher energy densities 
and temperatures (5Ml2|. 

In the literature various perturbative approaches are 
used for investigating (transport) properties of matter 
under extreme conditions [13--18J. In this work we con- 
sider an interacting (isospin symmetric) pion gas and fo- 
cus on its shear viscosity rj. Therewith we approach the 
properties of hot QCD matter from temperatures below 
the crossover region from hadronic to partonic matter: 
T < 140 MeV. A standard method to deal with such 
problems was developed by Kubo relating dissipative 
quantities to retarded correlators [TU] • We consider a hot 
pion gas close to thermal equilibrium, so we first perform 
an expansion in the dissipative forces d^F v . Then we im- 
plement a skeleton expansion of the four-point function 
entering the thermal viscous correlator Hp(ijj n ). The cor- 
relations between full propagators are neglected, which 
results in a quasi-particle approximation. Finally, the 
chiral expansion is used for calculating the pion self en- 
ergy at finite temperature. We find that the ratio r]/s 
(with s the entropy density) decreases as function of tem- 
perature but always exceeds the AdS/CFT bound 1/47T. 



II. RELATIVISTIC HYDRODYNAMICS 

We start with a macroscopic view on dissipative flu- 
ids. By definition we are dealing with a fluid that 
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does not fulfill Pascal's law, i.e. its pressure exerts 
transverse surface forces. Classically the dynamics of 
such a fluid is described by the Navier-Stokes equation. 
In relativistic hydrodynamics one extends the energy- 
momentum tensor of a perfect fluid by the dissipative ten- 
sor t^ v . With the Lorentz factor 1/7(2) = y/T— v 2 (x) 
let u^{x) = , y(x)(l,v(x)) denote the four-velocity, e(x) 
the energy density and P(x) the pressure. Then the 
energy-momentum tensor reads 



= u?u"(e + P) - PtT + t" 



(1) 



The form of the dissipative tensor t^ v can be fixed by 



the following three conditions [2Q] [2~T]: 



(a 



direct consequence of ([Tl when evaluated in the rest 
frame), <9 M s M > (second law of thermodynamics with 
the entropy-density current s M = su^\ and the assump- 
tion that only first-order derivatives d^v? contribute to 
t^ v ' . In this case the dissipative tensor is parametrized 
by the shear and bulk viscosities, r\ and £, respectively: 

t"" = ry^jV + fljX - -A""^ ■ u)] + CA^(9j_ • u) , 

(2) 

where we have introduced the following notations 

A" v = g" v - u^u u , d^ = d^- u»{u ■ d) = A^d v . (3) 

The second law of thermodynamics ensures that both 
transport coefficients are non-negative: r), £ > 0. We 
concentrate on the shear viscosity which describes the 
traceless part of t^ v ' . A more detailed discussion of the 
hydrodynamical properties with inclusion of shear and 
bulk viscosity can be found in Refs. |20[ 122] . 

Relativistic hydrodynamics is applicable to dissipative 
systems close to global thermal equilibrium. Experiments 
at RHIC suggest that the hot matter produced in heavy- 
ion collisions reaches local equilibrium after a short ther- 
malization time, t < lfm [TH71 123"]. A system in local 
equilibrium is characterized by the fact that one can di- 
vide the macroscopic system into mesoscopic zones in 
such a way that each zone by itself is in thermal equilib- 
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The more dissipative the system is, the more it devi- 
ates from global equilibrium: four-velocity, temperature, 
energy density, pressure, etc. display stronger space- 
time variations. The concept of local equilibrium there- 
fore translates into small values of gradients, e.g. d^u v . 
Therefore second-order terms can be neglected as already 
done in the parametrization of t» v in Eq. (j2j. This ap- 
proximation is justified when the system is still dissipa- 
tive but close to global equilibrium. The matter produced 
at RHIC meets this condition, after thcrmalization. 



equilibrium systems, po = p\^ =0 ■ Note that H is not a 
free Hamiltonian but represents a fully interacting theory. 
We mention other methods, e.g. the Boltzmann equa- 
tion, for investigating non-equilibrium systems 1 131 1 15] • 

The following quantities are introduced in order to de- 
compose a given energy-momentum tensor T pv : 



= A^ p u a T pa , x 



P=--A^T»\ 
— (A A - -A A )T pa 



III. NON-EQUILIBRIUM THERMODYNAMICS with ^ § iven in Ec l- ©> so that: 



The shear and bulk viscosities, r\ and f, are macro- 
scopic parameters for the dissipative part, Eq. ([2]), of the 
energy-momentum tensor. Now we recall the Kubo-type 
formula for the shear viscosity following Zubarev's ap- 
proach [2~4l |2"5] . First we introduce a statistical operator 
in the Schrodinger picture: 



1 

Q 



exp 



d 3 x B(x) 



(4) 



where Q = Tr exp [— J d 3 x B(x)] ensures Tr p = 1. We 
have defined the time-independent operator, 



B(x)=F»(x,t)T QlM {x,t)- 



dt'T lil/ (x,t')d>*F u (x,t'), 



with d t B(x) = 0, and the four-vector 

F"(x,t) =/3 B (x,tK(x,t). 



(5) 
(6) 



Here we have introduced the inverse proper temperature, 
ft = ~ = I , (7) 



again with the Lorentz factor 7. Evaluating /^.T 1 in the lo- 
cal rest frame (the distinguished frame of the heat bath) , 
one recovers the common temperature T . Note that with 
p s a Lorentz scalar, is indeed a four- vector. The oper- 
ator B is Lorentz variant, but with the additional spatial 
integration, the statistical operator p is indeed a Lorentz 
scalar. With f3 — l/T being the inverse temperature, it 
can be written as 



1 



exp[-£ff + A] 



(8) 



where f3H = J d 3 x F^(x,t)T ^(x,t) and A denotes the 
spatial integral over the second part of B(x) in Eq. 

The statistical operator p in Eq. ^ provides one pos- 
sible approach to non-equilibrium thermodynamics. It is 
not mandatory to introduce an operator B as done in Q , 
but the resulting form for p is physically meaningful: first 
because it depends on the dissipative forces d p F u , a mea- 
sure for deviations from global equilibrium. Secondly, 
for vanishing A, it reproduces the statistical operator for 



PA^ + K»u v + K v u^ 



nv 



(9) 



(10) 



Within linear response theory one can derive the impact 
of the dissipative forces on the energy-momentum tensor. 
We assume these forces to be small compared to typical 
energies of the system: (f3H) 3> (A)o, where (-) de- 
notes the thermal expectation value with respect to the 
equilibrium statistical operator po, i.e. (-)o = Tr(-po). 
Up to linear order in the dissipative forces the statistical 
operator in Eq. ^ reads 



Po 



1+ fd(e^Ue-^-(i) ( 
Jo 



(11) 



In particular, the linear response of the microscopic 
viscous-stress tensor ■K pv to dissipative forces leads to 
a connection between its correlation function and the 
macroscopic shear- viscosity parameter: 

V&t) = ^J dV J df (ir^(x,t),ir^(x" ,t')), 

(12) 

where the structure of the correlator (X, Y) follows from 
the statistical operator (111: 



{x,y) = ±J\t (x [e^y e -^-<y) ]> c 



(13) 



One can express this correlator as a real-time integral 
over a retarded correlator: 

(7r^(xU^(x',t')) w-i f dt(^{x),n^{x',t)) R , 

p J — 00 

(14) 

where x = (x, t) and the retarded correlator is defined as 
(X(t),Y(f)) R = -i0(i-f )<[*(*), F(f)]> • (15) 



The approximation (14) becomes exact in the large-time 
limit t — > 00, when the system reaches global equilib- 



rium. Finally, when combining Eqs. (12 I and (14), and 



evaluating at the local rest frame (/3 S /P — > 1), one finds 
the Kubo-type formula for the shear-viscosity field: 



1 



V (x,t) = / dV / df / dt (7T^(x,t), TV, (£',*)) 



R ■ 

(16) 
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Details concerning this derivation of r\ are worked out in 
Ref. [H]. For completeness we also give the expressions 
for the bulk viscosity £ and the heat conductivity k: 



C(x,t) = - / dV / dt' / dt(P'(x,t),P'{x',t)) R , 



CO J — CO 



(17) 



K (x,t) = - I dV / dt' / dt{K ll (2,t),K»(3! , ,i))n, 



CO J — CO 



(18) 

where P' = e d g^° —P- The heat conductivity can be ne- 
glected [26 if the chemical potentials are small compared 
to the temperature, i.e. if fJ-j/T -C 1. 



where is a real scalar field, m denotes the bare particle 
mass and g > is a small coupling constant. The viscous- 
stress tensor is determined entirely by the momentum- 
dependent parts of the Lagrangian: 



1 



(20) 



In general, the shear viscosity and other quantities such 
as pressure, temperature, four-velocity, etc. are fields, 
but on large time scales, when local equilibrium ap- 
proaches global equilibrium, they reduce to temperature- 
dependent numbers. We evaluate the shear viscosity ( 16 1 
at the origin of Minkowski space: 



7,(0) 



1 

To 



dt' / dtn R (t), (2i) 



IV. GENERAL PERTURB ATIVE EVALUATION 
OF THE SHEAR VISCOSITY 

A. Skeleton Expansion 

In this section we show how an expansion of the four- 
point correlator in full Matsubara propagators offers a 
suitable method for calculating the shear viscosity. The 
conccptional framework for this expansion is first worked 
out in gcj) 4 theory: 



C=\{d^){d^)- l -m 2 ci> 2 



(19) 



with the spatially-integrated retarded Green's function 
for the viscous-stress tensor 



II R (t) = -i J dV ([vr^(0),^(f',t)]) . (22) 

Using the analytical continuation to Minkowski space 
via iw„ >->• po + ie, we need to calculate the spatially- 
integrated thermal Green's function, labeled by /3: 



n /J (r)= / dV (T T [^(0)^(f',r)]) o 



(23) 



where 7~ T stands for the time-ordering prescription. Con- 
sider next the Fourier transform of Ilg (r): 



n^(w„) 



V 2 



dr e u 



d s p 



(amp a- _ A"W) Pli p v p p p a (T T [0(O)<KO W, rW-p, r)]) , (24) 



(2tt) 3 v ~ ~ 3 



V 2 



(T T [0(O)^(O)0(p,r)0(-p,r)]) o 



where <P(p,t) — jd 3 x e^' x 4>{x,T). The volume factors order, when neglecting these correlations, one obtains 
come from residual Fourier integrals originating from the 
normalization V J d 3 p = (27r) 3 . They ensure the correct n 
mass dimension of the momentum-integrated thermal 
Green's function: dim 11,3 = 4, since dim</>(p, r) = 
—2. The periodic boundary conditions of the bosonic 
fields in imaginary time are realized by the compact in- 
terval t £ [Q,0\ in Eq. pi} . 

Consider now the four-point correlation function 
(T T [<f>(0)<fi(0)<j)(p, r)<f)(—p, r)]) in detail. Let us recall the 



V 2 



(% [mm r)}) (% mm-p, r)]> 



(26) 



2 (Tr [mm)) (Tr \m r)H~P, T)]) 



covariance of two random variables X and Y: 



_((r T [0(O)0(p,r)]) o ) 2 =2&0,t) 



Cov(X,Y) = ((X-(X)) (Y-{Y)) ) = (XY)-(X)(Y) . 

(25) 

It demonstrates that the expectation value of a product 
factorizes if and only if the random variables are uncorre- 
cted. In this sense we introduce the skeleton expansion 
as an expansion in full propagators subject to correlations 
between them treated perturbatively |241 127| . At leading 



The first approximation comes from forming all possi- 
ble contractions that lead to full propagators (two-point 
functions) and factorizing them. The second one comes 
from neglecting the vacuum loops. Moreover the bosonic 
Matsubara propagator does not dependent on the direc- 
tion of the momentum flow: Gp{p,r) — Gg{— p, r). The 
skeleton expansion in gcj) 4 theory can be represented up 



4 



to next-to-leading order by the following diagrams: 



1 

V 2 



(r r [0 2 (o) 2 (p»]) c 




Here, the lines denote full bosonic Matsubara propaga- 
tors with self-energy insertions included. In the following 
calculation we restrict ourselves to the one-loop diagram 
since the primary application will be focused on pions 
within chiral perturbation theory. However, in scalar g0 4 
theory a resummation of ladder diagrams |281 |2"9"] is nec- 
essary to obtain the correct leading-order result for small 
coupling constants g -C 1. We briefly discuss this resum- 
mation and its impact on chiral perturbation theory in 
Appendix [B| 



Evaluating Eq. ( 27 1 at one-loop level leads to 



d 6 p 
(2tt) 3 



(28) 



where we have abbreviated the contraction of p 4 in 



It is convenient to express the Matsubara propagator 
in a spectral representation, 



Gp{p,r) 



du> p{ui,p) n(u>) e' 



,|r| 



(29) 



with the periodicity Gp{p,r) = Gp(p,T — 0). Here 
n(uS) 



where we have introduced the renormalized energy E p 
\Jp 2 + TOq + Re £r(p) and the spectral (half) width 



lip) = 2 r (P) 



2 P( 



Im£ R (p) 



(31) 



The spectral function p{p) = p(po,p) is determined by 
the imaginary part of the full retarded propagator, 



Pip) 



1 

2tti 



1 



1 



iPo - hip)f ~ El (po + ijipjf-El 



involving the difference between advanced and retarded 
propagators. We have used the identity —np(p) — 
Im G R (p). 



modified Breit-Wigncr form: 



The spectral function (|32j) can be written also in a 

(33) 



Pip) 



1 



* iPo - l 2 ip) - ED 2 + 4p2 7 2( p ) 



In this respect our treatment is equivalent to a quasi- 
particle approximation. 



B. Analytical Continuation 



l)" 1 denotes the Bose distribution and Substituting the Matsubara propagator (|29| into Q28) 



p(cu,p) the spectral function. 

On the other hand the retarded full propagator includ- 
ing the complex self-energy E R reads: 

G R » = p 2 - ml - ReS R (p) - iImE R (p) 

^( Po +ijip)) 2 -E 2 '" , ' 



the Green's function of the viscous-stress tensor can be 
expressed in a spectral representation as well. Using the 
temporal periodicity of Gp{p, t), the r-integration with 
a subsequent analytical continuation from discrete Mat- 
subara frequencies, cj n , to continuous energies, po, via 
iuj n t— > p Q + ie leads to 



n R (p ) 



d 3 P 
W) 



f P (P P V) P^Pv) 



dwi 



dui 2 p(oJi,p)p(u> 2 ,p) n{u\)n{w2) W e (wi2,£>o) , 



(34) 



where we have introduced the sum uj\2 = lo\ + W2, and This crucial property of W e guarantees a non-zero value 
the po-dependent factor W e : of the shear- viscosity (21 ): 

1 



1 



W t (Lu 12 , Po ) = — — — — . (35) 

Po + ie - con Pa + ie + w i2 

The real part of W e is symmetric in po, 

ReW,(- Po ) =ReW.( Po ) , (36) 



, ° dt' /' dt r ^ e-^ * n R ( Po ) 
, dt / dt' 



i 



2tt 



e~ ipo * n R ( Po ) (38) 



i 



10 d Po 



n R (p ) 



po=0 



i •-, , -, . , , n j where we have first interchanged the order of the inte 

while tne imaginary part changes its pole structure under . ,~ ,. ° . 



this reflection: 



gration (f, t') and then used the functional identity 



Im W e (-p ) = Im W-e(po) 



(37) 



r° d 

dt dt' e- ipot h> -2ni S(p ) 

-OO J t 



(39) 
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If Hr were an even function in po, the po-derivative eval- 
uated at po = would simply vanish. Since all the pq- 
dependence of ITr is carried by W e , its pole structure 
generates a non- vanishing, but highly singular expres- 
sion: 



lOdpo 



Po=0 



(40) 



Accordingly, only the diagonal line, W12 = u)i + 1^2 = 0, 
contributes to the shear viscosity. Denoting the perpen- 
dicular coordinate by u> — |(cJi — W2)) we obtain 



77 = 



57T 



d 3 p 



3 P^ L P V) P( P P V ) 



dio F(u) , (41) 



with the integrand: 



/?7 2 



(e#- - 1)" [£2 _ ( w _ i7 )2 



(w + i7) 2 



There are infinitely many poles of F(oj) in the complex 
plane (Fig. [TJ: the exponential function generates poles 
at uj = io; n , for n € Z \ {0}, where w„ = 27rinT are 
bosonic Matsubara frequencies. Indeed, the integrand 
-F(w) is regular at uj — 0. Furthermore, there are four 
poles (of second order): ujj = ±E p ± 17, for j = 1, 2, 3, 4. 
None of these poles is located on the real axis if 7 > 0. 
We use residue calculus to evaluate the integral, closing 
the contour in the upper half plane: 



/oo 
duj F(u) = 27ri 
-OO 



^ResF(iw„) + ^Res^(wj) 



.71=1 



(43) 

The two residua at W1.2 = ±-E p + i7 can be calculated 
analytically (with some lengthly result) and their con- 
tribution to the shear viscosity can be expanded in a 



00 



o 



o 



E r + i-y 



i7 o 



o 



-4mT 



-67riT 



Figure 1. Pole structure of the integrand F(oj) in Eq. (41 1 



For small 7 > the Sp-dependent poles are close to the real 
axis. The Matsubara poles on the imaginary axis give no 
contribution to the shear viscosity at leading order in 7. 



Laurent series in 7: 



27ri^ResF'(u; J 



7T/3 



( e /?£p - if 16£2 7 



+ 0(7). (44) 



Note that there is no constant term in the expansion. The 
spectral width Y = 27 is assumed to be sufficiently small 
for this expansion to hold. In view of the perturbative 
origin of Y this assumption is justified. 

Furthermore, one finds for the contribution of the n-th 
Matsubara frequency the following asymptotic behavior: 



Res F(iu) n ) 



n 9 



for small 7, f3 and large n . (45) 



The sum of ResF(iw„) over n G IN converges and yields 
a term of order 0(j 2 ) that can be neglected for small 7 
compared to the dominant 7 -1 contribution in (44 1. 



Evaluating the integral (41 1 in the local rest frame to 



leading order in Y, we arrive at the formula 



30tt 2 



dpp 



n{E p )[l + n{E p ) 



0(Y) . (46) 



This is the final result for the shear viscosity r\ in g<f) A 
theory. Note, that since Y(p) > 0, also rj > follows, as 
expected from the second law of thermodynamics. 

If one turns off the interaction and considers a free the- 
ory, the spectral width Y degenerates to a delta function, 
resulting in a divergent shear viscosity. This seemingly 
counter-intuitive behavior has been known for a long time 
|30| . not only for bosonic systems |31) . 



V. SHEAR VISCOSITY OF A HOT 
INTERACTING PION GAS 

A. General Discussion 

Next we apply the skeleton expansion to pions within 
the framework of chiral perturbation theory (%PT). For 
T < 160 MeV [32 36 , i.e. in the temperature range where 
chiral symmetry is spontaneously broken, X-PT can be 
applied with confidence. We use the second-order chiral 
Lagrangian, £2 (in c-gauge) , expanded up to fourth order 
in the pion fields: 



+ i(*-0/x*)(*-0"*) 



72 I 77 ) + 



(47) 



In contrast to the discussion in g0 4 theory, not only the 
kinetic term, but also the momentum-dependent interac- 



tion terms contribute to the viscous-stress tensor ( 20 ) 



1 



7i — ('A P A <T A A pcr 



(d p TT ) • (<3 CT 7T ) + ^— (7? • <9 p 7T ) (7T ■ <9 CT 7? ) 

J 7T 



(48) 
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The momentum-integrated thermal Green's function in 
Fourier space now has additional terms in comparison to 



g<f> A theory. The kinetic term is analogous to Eqs. (24) 
and f26|: 



A = ^2 ^ ' ^(°) r ) ' r )l)o ■ ( 49 ) 

In addition there exist two terms at order 0(f~ 2 ), and 
three terms at order 0(f~ 4 ): 



B l = 



Bo = 



1 



f%v 4 



(7;[#(0)7?(-^T)-7r(?,r)]> ( 



(50) 



Green's function Hp(u) n ) including corrections of order 
C(/~ 4 ). In the thermodynamic limit. V — > oo, this cor- 
rection term vanishes. Thus we find for the shear vis- 
cosity in %PT the same functional formula rj[T] as in g(f> A 
theory, Eq. (46 1, but corrected by an isospin factor three: 



P 



IOtt 2 



dpp 



6 n(E p )[l + n(E p )} 
E*T{p) 



(52) 



with E 2 — p 2 + m 2 . We mention that this identity does 
also hold if one allows for higher-order terms in the chiral 
expansion or higher-order terms in the field expansion of 
£2 ■ All these corrections would be of higher order in the 
coupling 1//^ and suppressed by higher powers of the 
inverse volume. 



B. Analytic Evaluation of the Spectral Width 



c 



<r r [^ 4 (0)7f 2 (-p,r) 7 r 2 (p ! r)]) 



The disconnected parts of the skeleton expansion are 
vacuum loops and can therefore be dropped. Note 
again that the diagrammatic representation already uses 
Gp(—p,T) = Gp(p,r). Furthermore, the first terms of 
Bi , B 2 and the second one of C can be absorbed by the 
skeleton vertex (or eliminated by a field renormalization) , 
since we are anyhow expanding in fully- dressed quanti- 
ties. At order C(/~ 4 ) there is an additional term in 
the skeleton expansion giving rise to a correction of the 
momentum-integrated thermal Green's function: 



TLp{u n ) = 2 (N 2 - 1) dr e 1 ^ J 



d 3 p 
(2tt) 3 
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(51) 



Additonally to the calculation in g<fi theory ( 28 ) we have 
now the flavor prefactor N 2 — 1 = 3 which arises from the 
sum over the isospin degree of freedom. The leading order 
term in Eq. (51 ) has the symmetry factor 2, whereas the 



correction term has the combinatorial factor 4! = 24 . 

We have started with the second-order chiral La- 
grangian and arrived at the skeleton-expanded thermal 



Once the spectral width T is known, the shear viscosity 
j] in Eq. ( [52] ) can be computed. In xPT the leading- 
order contribution to T comes from the two-loop diagram 
shown in Fig. [2] The self-energy diagram at one-loop 
level is real, hence it gives no contribution to the spectral 
width. 

In the two-loop diagram in Fig. [2] the two Matsubara 
sums over (mi, ra-i) £ Z 2 , can be performed analytically. 
Taking into account the energy-momentum conservation, 
E 2 = (p — pi — P2) 2 + m 2 , and the on-shell identity for 
the Mandelstam variables, s + t + u = 4m 2 , we arrive 
at an expression for the thermal self-energy Y,p{iuj n ,p) 
which can be analytically continued via ico n 1— > po + ie 
to the retarded self-energy Sr. The spectral width of a 
pion with four- momentum p^ = (po,p) reads: 




Figure 2. Self-energy diagram of the pion in the heat 
bath. The arrows indicate the direction of the three- moment a. 



r(p) 




4/4 n ( Ep )E, 



2(s 2 +t 2 



-u 2 )-9mi]n(E 1 )n(E 2 )[l + n(E 3 )]{27r) 4 S^(p- Pl -p 2 +p 3 ) . (53) 



The combination of relative signs in the delta function 
§( 4 )(p — p 1 — P2+P3) is the only one that gives rise to an 
on-shell contribution to the imaginary part Im Sr. In all 



other cases the energy-conservation delta function leads 
to zero. In general, i.e. in the off-shell case, there are four 
different combinations of relative signs, shown in Fig. [3] 
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The first case leads to 5{pv+ E\+ E2+ E3) = 0, hence this 
realization gives no contribution. The structures ( — h+) 

and ( ) are not realized on-shell, because they would 

describe the decay of a particle into three copies of it. 
This is not possible for massive pions. Thus, only three 

equal terms with ( h) signature remain. 

From a numerical point of view the representation of 
r in Eq. (53 1 is only symbolic. The (four-dimensional) 



delta function and the nine-dimensional integral over it 
are not well-conditioned for numerical evaluation. As 
shown in the appendix [A] one can rewrite the expression 
for the spectral width as a two-dimensional energy in- 
tegral suitable for numerical evaluation. This procedure 
introduces the quantities (here pi = \pi\) 



Q_ = maxj>2 — P3I, \P-Pi\} 
Q + = min{p 2 + P3,P + Pi} ■ 



(54) 



The integrand in Eq. (|56| depends on differences of pow- 
ers of Q± , 



X(Q-,Q, 



E 

ie{0,±2,±4} 



(55) 



where the five coefficient functions Xi have the mass di- 
mension dim Xi = 4 — i and depend on the pion mass 
and on the energies E2 , E3 , E p on ly. Their explicit forms 
are given in the appendix, (Al8l to ( A21 1. In the end 



the double-integral representation of the spectral width 
reads 



r(p) 



sinh 



f3E p 



dE* 



2pE p f* (8tt)3 

X(Q +: Q_)e(Q+ 



E p +E 3 -r 



dE 9 



(56) 



sinh^fi 



sinh 



2 



sinh 



HE* 
2 



In comparison with Ref. |37| . where the mean free path 
of hot pions has been calculated in the infinite-volume 
limit (consistent with our skeleton-expansion approach), 



-<- 



(+ + +) 



(- + +) 



->- 



(-- 



(- 



Figure 3. Relative signs of energies in the two-loop diagram 
with respect to the incoming-particle energy po. Only the 
2 — > 2 scattering process with signature ( (-) gives an on- 
shell contribution to the spectral width T(p) in Eq. (531. The 



dashed lines symbolize the on-shell condition at the vertices. 
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'*••.. T=93MeV 




T = 108 MeV 






/ 
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■ / ^ 




:/ / 




11/ 






Figure 4. Mean free path A(p) of a hot pions for differ- 
ent temperatures (upper figure) and pion masses (lower fig- 
ure). The two solid lines in the figures are drawn for the 
same combination of mass and temperature: (m-n- , T) = 
(140 MeV, 124 MeV). 



we have found an analytical representation of the spec- 
tral width with a different algebraic structure. Our re- 
sult involves the rational function X(Q + , Q_) whereas a 
polynomial function of the energies appears in Ref. |37| . 
In addition, the max/min pattern in the two represen- 
tations are quite different and it is not apparent how to 
convert one into another. 

Irrespective of that our results for the spectral width 
are identical with those of Ref. [37] : only the integrands 
differ, but the numerical values of the integrals for the 
mean free path coincide. In essence, the mean free path 
X(p) is the inverse spectral width: 



X(p) = 



P 



E p T(p) 



r(p) 



(57) 



As an aside, we mention that the spectral width T(p) in 
g<fi theory can be obtained by inserting the expression 
X = 576 f%g 2 (Q + — Q-) into Eq. (56 1 and replacing 



by m. This result for X is derived from the two-loop 
diagram in Fig. [2] as well, but evaluated in gcj> 4 theory. 
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Figure 5. Shear viscosity r](T) of an interacting pion gas for 
different pion masses. In the chiral limit one gets r\ ~ 1/T as 
expected from dimensional analysis. 



Figure 6. Entropy density s(T) of an interacting pion gas 
for different pion masses. We compare the xPT results to 
the Stefan-Boltzmann limit s(T) = ^j-T 3 of three massless 
non-interacting pions. 



C. Numerical Results for A and rj 

The numerical results for the mean free path A(p) 
are shown for different temperatures and pion masses in 
Fig. |4j For increasing temperatures the abundance of 
thermal pions lowers the mean free path, but the max- 
imum position p w 100 MeV is almost independent of 
the temperature. The mass dependence of X(p) is more 
pronounced: when switching from the physical value, 
TOjr = 140 MeV, to pions with half the physical mass and 
to the chiral limit, the mean free path features a mono- 
tonic increase at low momenta p < 200 MeV. For all 
temperatures and masses the mean free path decreases 
for large momenta. In the numerical evaluation we have 
used the value f % = 93 MeV for the pion decay constant. 

Interpolating the numerical results for the mean free 
path at different temperatures and masses finally allows 
to calculate the shear viscosity. The results are shown 
in Fig. [5j In the chiral limit, = 0, one finds a re- 
ciprocal dependence on the temperature, r)(T) ~ f^/T, 
as expected from dimensional analysis. For the phys- 
ical pion mass, m v — 140 MeV, the maximum of r\ is 
located at T « 20 MeV. With decreasing pion masses, 
this maximum moves to lower temperatures. For high 
temperatures the shear viscosity depends only weakly on 
the pion mass. 

Note, that the function ?7(T;m w ) is not continuous at 
the origin, since 

= lim lim r](T;m 7r ) ^ lim lim ^(T^rrO = oo . 

m„— s-0 T— s-0 T— s-0 m„—>0 

v ' > v ' 

= r^l/T 

n (58) 
In Fig. we have chosen the units MeV/fm 2 instead 

of MeVMn order to meet the classical interpretation of 

shear viscosity. 



VI. RATIO 7] Is FOR THE PION GAS 

In 1998 Maldacena (38] demonstrated that, under cer- 
tain conditions, there is a duality between superstring 
theory and superconformal field theory. More precisely, 
the AdS/CFT correspondence between supergravity on 
five-dimensional anti-de Sitter space, AdSs, and four- 
dimensional M = 4 superconformal SU(V r c ) Yang-Mills 
theory in the 't Hooft limit has been proven. In this 
special case the ratio rj/s (shear viscosity over entropy 
density) is equal to the KSS lower bound l/4vr [39]. In 
2005 the KSS conjecture has been formulated: 

Most quantum field theories do not have sim- 
ple gravity duals. (...) We speculate that the 
ratio rj/s has a lower bound rj/s > 1/47T for 
all relativistic quantum field theories at finite 
temperature and zero chemical potential. 

So far, the KSS lower bound for rj/s is respected by exper- 
imental results for a wide variety of thermal systems |40| . 
The fundamental theory describing heavy-ion collisions is 
QCD which does not possess a gravity dual: QCD is nei- 
ther supersymmetric nor conformal. The classical scale 
independence of the QCD Lagrangian is broken anoma- 
lously by quantum effects resulting in a running coupling 
a s (q 2 ). In addition, QCD is a SU(3) gauge theory and 
can be described by the large-iV c limit only approxima- 
tively. Since 2007 counterexamples to the KSS conjecture 
have been established (4"Tl 1412] . 

The temperature-dependent entropy density of the in- 
teracting pion gas reads at two-loop order |33| 

s(T) - [4T 2 h 5 (pm n ) + 3mlh 3 (pm n )] 

3?77 2 T 

- 77ri72 HPm w ) [2T 2 h 3 (pm w ) + m£/n(0m*) , 
16tt 4 /^ 

(59) 
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Figure 7. Ratio between shear viscosity and entropy den- 
sity ri(T)/s(T) of an interacting pion gas for different pion 
masses. The horizontal dashed line is the lower bound 1/4-7T 
from AdS/CFT correspondence. 



with the integral function for £ > : 



M0 



d.r 



(60) 



In the limit / w — > oo only the first line of Eq. ( 59 I con 



tributes to the entropy density. The correction to s at 
order 0(f~ 2 ) is less than 1.5% (for m % = 140 MeV) and 
decreases even more for lower masses. In the chiral limit 
the entropy density is given by 



lim s(T) 

n x -s.0 



27T 2 

15 



(61) 



This is the well-known Stefan-Boltzmann limit for three 
non-interacting massless pions. In Fig. [6] we compare the 
reduced entropy density s/T 3 of the pion gas in x?T to 
the constant 27r 2 /15 from Eq. ( |6l) . 

In Fig. [7] we show the temperature dependence of the 
ratio rj/s for the interacting pion gas. Actually, for tem- 
peratures up to T ss 140 MeV we find a ratio rj/s which 
is well above the KSS lower bound 1/4-7T. In the consid- 
ered temperature range 80 MeV < T < 140 MeV the ratio 
rj/s is almost independent of the pion mass. Our results 
compare qualitatively well with alternative approaches to 
calculate the shear viscosity of the pion gas [431 HI]. 



VII. SUMMARY 



In this work we have calculated, on the basis of the 



Kubo-type formula (16 1, the shear viscosity 77 of a pion 



gas within chiral perturbation theory (%PT). The skele- 
ton expansion, corresponding to an expansion in full 
propagators, has been applied to the thermal (four-point) 
Green's function. At one-loop level this expansion leads 



to the squared Matsubara propagator (two-point func- 
tion). This technique has first been explored in g<j) the- 
ory, and then applied to %PT. We have found that at 
one-loop level of the skeleton expansion and in the ther- 
modynamic limit V — > 00, the shear viscosities for real 
scalar fields, Eq. (46), and pseudoscalar pions, Eq. (52 I, 



differ only by an isospin factor. The functional depen- 
dence of rj on the spectral width coincides for both cases. 

The spectral width T(p) of pions has been calcu- 
lated to two-loop order in thermal xPT resulting in a 
new analytical representation Eq. (56). We have found 



that the shear viscosity of the pion gas reaches a max- 
imum at T « 20 MeV in the case of physical pion mass 
TOtt = 140 MeV. Furthermore we have investigated the 
ratio rj/s (with s the entropy density) and found that 
it decreases monotonously with rising temperature but 
exceeds the KSS lower bound 1/4-7T for all temperatures 
T < 140 MeV where xPT is applicable. 



Appendix A: Details on the Analytical Evaluation of 
the Spectral Width 

We outline the analytical calculation which leads to 
the numerically well-conditioned expression ( [56]) for the 
spectral width T(p). The momentum integral J d 3 p\ in 
Eq. ( 53 1 is canceled by the three-dimensional delta func- 



tion, hence we can express the corresponding energy E\ 
by the incoming momentum and the two remaining in- 
ternal momenta: 



El = E 2 + \p 3 - p 2 \ 2 + 2p\p 3 - p 2 \ cos 9 P 



(Al) 



We integrate out the dependences on cos 9 p and <p p , using 



1 = 



1 

2tt 



2tt 



dtfp 



1 



d cos 9 



p ■ 



(A2) 



Denoting the squared on-shell pion-pion scattering ampli- 
tude by V 2 n = [2(s 2 + 1 2 + u 2 ) — 9m*]/ f% and performing 
the cos#p integral, the energy-conserving delta function 
fixes the value of cos 9 n : 



d cos I 



V* n (cos (p p , cos 9 P ) 



2E 1 



5(E p ~E 1 -E 2 +E a ) 



ApQ 



(A3) 

where we have introduced the abbreviations Q = \p3 — 
p 2 \ and x = [(E3 - E 2 )(E 3 - E 2 + 2E p ) - Q 2 ] . In 
addition to 9 P = <(p,P3 — p 2 ) and ip p , we introduce the 
following quantities: 

(f = <{p,P3+p 2 ) , 

A = <(p 3 -P2,P3+P2), (A4) 
R = \P3 +p 2 | ■ 

It follows immediately that R can be expressed in terms 
of Q: 



R 



2E 2 



2E 2 



4ml - Q 2 



(A5) 
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Furthermore, it is possible to express V 2 n just in terms 
of the sum and the difference of p 2 and p 3 , using s — 
2m 2 ~ 2 b ■ P3)> t — 2m 2 ~ 2(p ■ Pi) and u = Ami - s-t: 

/ 4 V 2 n = lm\ + 16ml [E p {E 3 - E 2 ) - pQ cos P ] 

+ 12[E P (E 3 - E 2 ) -pQcosQpf (A6) 
+ 4 [E p (E 3 + E 2 ) - pR cos . 

The second argument of V^ 2 n (cos <p p , cos 9 P ), is already 
fixed by (A3): cosO p = Xq. The first argument, cos<p p , 



can be related to cos <p, by spherical trigonometry: 



cos ip — cos 9 P cos A + cos ip p Jl — cos 2 p y/l — cos 2 A . 

(A7) 



Using its definition in ( A4 I and the relation A A5| , cos A 



is a function of the energies E 2 , E 3 and Q only: 



cos A 



Ps, -P2 P3+ P2 E 2 - El 



R 



QR(Q) 



(A8) 



nomial in cos 

-2 



Inserting (A8) into (A6l, VA. becomes a quadratic poly- 



V on (cos ip p ) = A + A 1 cos ip p + A 2 cos ip p , (A9) 

with some complicated coefficient functions Ai depending 
on E p , E 2 ,, E 3 , run, and Q. Now, we can easily carry out 
the integration over the azimuthal angle tp p from ( A2 ) : 



1 f 1 

— / difp (A + A x cos ifp + A 2 cos 2 <p p ) = A + -A 2 . 

Z7T Jq Z 

(A10) 

Consider next the two remaining momentum integrals 
J d 3 p 2 and J d 3 p 3 coming from Eq. (53). We have to 



ensure that the fixed value of cos6> p lies in the allowed 
range | cos# p | = |aco] < 1- This constraint is implemented 
by the factor Q(F(Q)), with 

F(Q) = Ap 2 Q 2 - [{E 3 - E 2 )(E 3 - E 2 + 2E p ) - Q 2 f , 

(All) 

because \xq\ < 1 means 

_ l< (E 3 -E 2 )(E 3 -E 2 + 2E p )-Q 2 
2pQ 

which is equivalent to F(Q) > 0. As a function of Q 2 , 
F(Q) is just a concave-down parabola with F(Q) > 
for Q\ < Q 2 < Ql with roots Q\ = {p - pi) 2 and Q\ = 
(p + pi) 2 . Since Q > by its definition, we arrive at 
Qi = \p — pi\ and Q 2 — p+pi, the only two non-negative 
roots of F(Q). At this stage, we arrive at the expression 



d 3 p 2 fd 3 p 3 



E-2 



E, 



^ roc rEp+Es — m^ 

- dE 3 dE 2 



J dn 2 J dn 3 (a + 1 -a 2 ^j Q(F(Q)) . 

(A 13) 

The two remaining solid-angle integrals can be reduced 
to an integral over Q, using the relation 



Q 2 =pI+pI-2P2P3C0S< 



'23 



(A14) 



Inspecting the coefficients Aq and A 2 , one observes that 
they consist of a few positive and negative powers of Q 2 : 



' f dQ 2 f &Vl 3 P -^(av+ 1 -A 2 \q{F(Q)) 



4p 



Q 



8tt 2 r p2+P3 
4p 

2ir 2 



P 



/ dQ [5X 4 Q 4 + 3X 2 Q 2 + X - X_ 2 Q~ 2 - 3X_ 4 Q~ 4 ] 6(F(Q)) 

J\P2-P3\ 

Xi {Q% - Ql) + X 2 (Ql - Q 3 _) + X (Q+ - Q-) + X- 



(A15) 



Q+ Q- 



X. 



1 Qi 



1 

Qi 



Here the characteristic form of X in Eq. ( A17 1 with dif- 



ferences of powers of Q± has emerged. The step function 
Q(F(Q)) gives rise to a max/min pattern of the proper 
boundaries Q±, which read 



Q_ =max{|p 2 -P3\,\p-Pi\} , 
Q + = min{p 2 +p 3 ,p + p 1 } . 



We combine the five coefficient functions 
Xi(Ep,E 2 ,E 3 ]m v ) in the last line of Eq. ( A15 1 to 
a new function X(Q+, Q-): 

X(Q_,Q+)= ]T ^ (Q+ +1 - Q- +1 ) • (A17) 

ie{0,±2,±4} 



(A16) The explicit expressions for Xj are the following: 
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X, = 



10' 



X 2 = - [6ml -E 2 p - AE 2 2 - 4£f + E p E 2 - E p E 3 + 7E 2 E 3 ] 



(A18) 



X n = 15m? 



8ml (3E 2 E 3 + E p E 2 - E p E 3 - E 2 p - 2E% - 2E 2 ) 



2E p (E p ~E 2 + E 3 ) {hE\ + 2E 2 E 3 + 5£f) + (E 2 - E 3 f (7E 2 - 4E 2 E 3 + 7E% 



(A19) 



X_ 2 =4ml (E 2 - E 3 f (2E p E 2 



2E 2 p - E 2 2 - E 2 



2E P E 3 ) 



■ 2E p {E 2 - E 3 f [E p -E 2 + Ea) (7E$ + 10E 2 E 3 + 7E 2 ) + 2(E 2 - E 3 f (2E 2 + 3E 2 E 3 + 2E 2 ) , 



(A20) 



X, 



1 



(E 2 - E 3 f {2E p - E 2 + E 3 f {E 2 



E 3 



(A21) 



Putting all pieces together, we arrive at a well- 
conditioned double-integral representation for the spec- 
tral width: 

2n~ 1 (E ) r°° rE p +E3-m* 

m = ^Jz p L I dE 3 l dE 2 



pEpf* (8tt)3 

• X(Q + , Q_)n{E x )n{E 2 )[l + n(E 3 )} 6(Q+ - Q_) . 

(A22) 

Alternatively, by using the identity 



ie-i^ sinh- 1 ( pEi 



(A23) 



and energy conservation E\ = E p — E 2 + E 3 , the spectral 
width can be expressed in the form given in Eq. ( 56 I . 
The additional factor 0(Q+ — Q_) in Eq. ( A22 1 must be 



introduced in order to restrict the double integral to the 
physical kinematical region in the i^-E^-plane where the 
upper limit Q + is actually larger than the lower limit Q_ . 
At the same time this guarantees that the spectral width 
r(p) is made up from strictly-positive contributions only. 



Appendix B: Discussion of ladder-diagram 
resummation for the skeleton expansion in xPT 

In order to approximate the four-point correlation 
function that enters Hp(u) n ) in Eq. (24) we have per- 
formed a skeleton expansion ( |27| from which we have 
taken only the one-loop term into account. A detailed 



analysis of this expansion for a real scalar theory with cu- 
bic and quartic self interactions has been given in [2"8ll29| . 
It is concluded there that for a consistent treatment of the 
shear viscosity r\ at leading order in the small coupling 
constant g <C 1, a resummation of all ladder diagrams 
needs to be performed: 

(...) this means that higher loop dia- 
grams can be just as important as the one- 
loop contribution if they are sufficiently in- 
frared sensitive. 129 



This argument is based on the observation that in the 
limit of vanishing thermal width, r — > 0, the occurrence 
of pinched poles spoils the usual perturbative counting 
in powers of the small coupling g (compare Fig. [T]). The 
one-loop diagram in Eq. (27 1 scales as 1/r, whereas the 
two- loop diagram scales as g 2 /T 2 , and the n-loop ladder 
diagram scales as g 2n /T n+1 . Hence, there is a dimen- 
sioned scaling factor g 2 /V for every additional rung in 
the ladder-diagram expansion. Since the spectral width 
is T ~ g 2 , all ladder diagrams are of the same order 
0(g~ 2 ), and therefore need to be resummed. According 
to Refs. |45l 135] this can be done in an efficient way by 
using a two-particle irreducible effective action. 



Let us discuss to which extent and scalar gcf> A 

theory differ in this respect. In the chiral limit, m n = 

0, the pion-pion interaction is purely of derivative type, 

1. e. proportional to p 2 . In such a situation, the infrared 
singular l/T(p) terms resulting from the nearly pinching 
poles are compensated by momentum-dependent factors 
in the numerator. Inspecting the chiral Lagrangian (47 1 
one can identify the dimensionless coupling g =p 2 /2f£. 



The additional factor appearing at each higher order in 
the ladder-diagram expansion is g 2 /T(p) = p 4 / (4/*r(p)), 
which vanishes in the infrared limit, p — > 0. However, the 
additional chiral-symmetry breaking mass term in the 
chiral Lagrangian gives rise to a pion-pion interaction 
analogous to the vertex in g0 4 theory, with g = ml /8f 2 . 
Following the scaling arguments of |28| this feature may 
require the resummation of ladder diagrams. 

In general, xPT is far less infrared sensitive than other 
bosonic field theories. While a detailed analysis needs yet 
to be performed, one expects that the numerical conse- 
quences of such a resummation may be less important 
than in gcf) 4 theory. Let us finally note that xPT is only 
applicable for low temperatures T < m^ = 140 MeV. 
In that temperature range, thermal corrections to the 
pion mass are smaller than 5% and therefore negligible 
[471 148| . For instance, this is manifest in the absence of 
a Linde problem [IS] in xPT. 
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